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INTRODUCTION 


In  many  situations,  the  operating  environment  of  an  active  pulsed  system 
(radar,  communications,  or  sonar)  limits  the  number  of  pulses  transmitted. 
Moreover  the  pulses  transmitted  through  a  medium,  reflected,  and  returned  can 
be  mathematically  modeled  as  received  pulses  with  pulse-to-pulse  amplitude 
variation  (fading).  An  air  traffic  control  radar  In  a  mountainous  region 
operating  during  severe  weather  might  be  subject  to  such  limitations.  High 
frequency  radio  communication  at  long  ranges  (skip  distances),  when 
Ionospheric  conditions  are  changing  because  of  normal  diurnal  processes  or 
during  solar  storms,  provide  a  second  class  of  examples.  Replacement  of  these 
electromagnetic  systems  by  their  acoustic  analogues  provides  the  further 
Instance  of  a  sonar  operating  among  Islands  or  at  convergence  zone  ranges 
(underwater  skip  distances).  The  signal-to-noise  ratio  at  the  system  input 
could  be  randomly  reduced  to  a  value  so  low  as  to  preclude  acceptable 
detection  performance  from  processing  each  signal  pulse  independently. 

Detection  performance  of  a  processing  system  can  be  mathematically 
characterized  by  the  parameters  Pg,  probabllity-of-detection,  and  Pp, 
probabi 1 Ity-of-false-alarm.  When  several  pulses  are  transmitted  in  a  pulse 
train  burst,  and  the  returned  echoes  are  processed  as  a  group,  we  wish  to 
determine  the  minimum  required  input  signal-to-noise  ratio,  as  a  function  of 
the  number  of  pulses,  K,  in  the  burst,  necessary  to  achieve  the  preassigned 
acceptable  level  of  performance,  pD»pp-  Similarly,  if  only  a  fixed  total 


energy  Is  available  for  transmission,  we  want  to  select  the  "best" 
distribution  of  the  available  energy  Into  K  equal  duration  pulses  to  ensure 
acceptable  system  performance.  Stated  differently,  what  Is  the  optimal  number 
of  baskets  (pulses)  for  our  eggs  (energy)? 

For  our  purposes  here,  a  transmitted  burst  consists  of  K  identical, 
nonoverlapping  (in  time),  sinusoidal  pulses.  The  pulse  train  received  is 
mathematically  modeled  as  time-delayed,  randomly  phase-shifted,  versions  of 
the  transmitted  pulses  with  random  pulse-to-pulse  amplitude  variations  of 
known  first-order  probability  density.  This  received  signal  is  further 
corrupted  by  additive  Gaussian  white  noise.  The  fading  is  assumed  to  be 
unnoticeably  slow  during  a  pulse,  so  that  the  only  fades  it  is  necessary  to 
model  are  from  pulse-to-pulse,  which  can  be  correlated  to  an  arbitrary  degree. 

One  possible  detection  tactic  might  be  to  transmit  the  smallest  energy 
signal  just  necessary  to  achieve  the  required  average  detection  performance, 
and  accomplish  it  with  a  single  pulse.  Because  of  a  variety  of  limitations, 
either  physical  or  financial,  this  tactic  may  not  achieve  the  desired 
performance.  Multipulse  bursts  might  work  in  cases  where  the  available 
slgnal-to-nolse  ratio  is  low.  Here  we  consider  in  what  sense  this  intuition 
Is  correct  for  one  class  of  correlated  fading,  and  we  compare  these  results 
with  system  performance  in  a  nonfading  medium. 

In  the  next  section,  the  characteristics  of  the  transmitted  and  received 
signals,  the  corrupting  additive  noise,  and  the  receiver  processing  are 


mathematically  modeled.  This  provides  notation  and  scenario  for  the  analysis 
of  the  hypothesized  receiver,  depicted  In  figure  1.  The  class  of  power  scaling 


Figure  1.  Receiver  Block  Diagram 

fading  variates  considered  have  chi-square  densities  so.  In  a  sense,  the 
received  envelope  amplitude  scalings  can  be  thought  of  as  the  square  roots  of 
chi-square  variates.  The  possibility  that  several  pulses  in  a  pulse  train  are 
correlated  Is  covered  as  well.  Numerical  examples  for  several  combinations  of 
P0  and  Pp,  for  some  representative  received  scaling  probability  density 
functions  and  pulse-to-pulse  correlations,  are  presented  and  discussed. 

Two  appendices  provide  mathematical  details  supporting  the  main  analysis 
by  deriving  the  characteristic  function  and  exceedance  distribution  of  the 
system  output;  a  third  appendix  lists  computer  programs  for  the  computation  of 
required  input  signal-to-noise  ratios  for  fading  and  nonfading  received 
signals . 
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PROBLEM  DEFINITION 


Signal  Characteristics 

The  detection  problem  of  Interest  Is  described  In  this  section.  A 
high-frequency  pulse  of  length  L  seconds  and  carrier  f  Hz  Is  utilized  as 
the  basic  signal  waveform  component: 

p(t)  «  cos(2*fQt  +  ®Q)  for  0  <  t  <  L  .  (1) 

We  have  LfQ  »  1;  that  Is,  each  pulse  contains  many  cycles  of  the  carrier. 
The  actual  transmitted  signal  waveform  Is 

K 

A^  p(t  -  tk)  ,  (2) 

k-1 

where  A  Is  the  transmitted  signal  peak  amplitude,  K  Is  the  number  of  pulses, 
and  {tk}  are  the  transmission  times  which  are  arbitrary  except  that 

tk+l  ~  >  1  ;  i,e-*  no  overlaP  •  (3) 

The  received  signal  Is 

K 

a2  rk  p(t  -  tk  -  D)  , 


(4) 


where  {r^  are  amplitude  scalings  Imposed  by  propagation  and  attenuation 
through  the  medium.  Bulk  delay  0  Is  common  to  all  pulses,  and  p  denotes  a 
randomly  phase-shifted  version  of  (1)  by  Independent  uniformly  distributed 
phase  shifts  {©k}  for  1  <  k  <  K. 

The  average  received  signal  power  during  a  single  pulse  Is 

A2  r2  cos2(2*fQ(t  -  tk  -  D)  +  eQ  +  0k)  *  1  A2  r2  .  (5) 

and  the  average  received  signal  energy  per  pulse  Is  then 

\  A2  r2  L  .  (6) 

This  Is  a  very  Important  parameter  regarding  the  received  signal;  It  Is 

-7  T 

Independent  of  pulse  number  k  under  the  assumption  that  rk  *  r  ,  that 
Is,  stationary  fading. 

Noise  Characteristics 

* 

There  Is  also  assumed  present  at  the  receiver  Input,  additive  zero-mean 
white  Gaussian  noise  n(t)  with  a  double-sided  spectrum  level  watts/Hz,  or 
equivalently,  a  single-sided  spectrum  level  Nq  =  2Nd.  The  total  received 
waveform  Is  therefore, 


Receiver  Processing 


The  bulk  delay  0,  as  well  as  the  time  separations  between  pulses, 
are  assumed  known  at  the  receiver;  otherwise,  a  search  over  time  delay  D  Is 
required.  In  addition  to  a  possible  Ooppler  search.  In  general.  The  receiver 
consists  of  K  matched  filters,  one  for  each  of  the  signal  pulses  In  (4);  these 
filter  outputs  are  envelope  detected  and  sampled  at  the  appropriate  Instants 
corresponding  to  the  peak  signal  outputs.  Finally,  these  envelope  samples  are 
squared  and  summed  over  the  total  of  K  pulses.  This  sum  Is  then  compared  with 
a  threshold  for  decision  about  signal  presence  or  absence;  see  figure  1. 

The  analytic  problem  of  Interest  Is  to  determine  the  probability  that  the 

decision  variable  above  exceeds  a  threshold.  From  this  quantity,  we  can 

determine  the  detection  and  false  alarm  probabilities  as  functions  of  all  the 

various  parameters  presented  In  (1)  -  (7).  Then  for  specified  detection  and 

false  alarm  probabilities,  the  required  Input  slgnal-to-noise  ratio, 

ET/N  ,  can  be  determined  numerically, 
i  o 
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ANALYTIC  RESULTS 

Characteristic  Function  of  Decision  Variable 


The  derivations  of  the  statistics  of  the  decision  variable  y  are 
conducted  In  appendices  A  and  B.  Specifically,  the  characteristic  function  of 
y  Is  given  in  (A-31)  as 

(l  -  152o^  ”  ,  (8) 


where 


i  V  • 


<*‘(1  +  R)  . 


e  -1  £ 
Nq  N  * 


and 


m  average  received  signal  energy  per  pulse _ 

Nq  *  single-sided  received  noise  spectral  density  level  •  *  ' 


N  -  rn  Ke,  Ke - £ -  .  (10) 

k.Jl-l 


Here,  m  Is  a  parameter  In  the  first  order  probability  density  function  of 
amplitude  scalings  {r^j ,  namely, 

p  ( u)  =  - exp( -u  /a)  f0r  u  >  0  (m  >  0)  .  (11) 
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To  explain  the  other  parameters  In  (10),  we  use  the  power  scaling  variable 


flk  ■  rk 


(12) 


for  the  k-th  pulse  In  (4)  and  (7).  From  (11),  the  probability  density 
function  of  {q^  Is 


,m-l 


p„(u)  -  u  exP(-fl/g)  for  u  >  0  (m  >  0)  . 


T(m) 


(13) 


This  Is  recognized  as  a  chi-squared  variate  of  2m  degrees  of  freedom;  however, 
2m  need  not  be  an  Integer.  It  Is  readily  verified  that 


Std  Oev(q) 
Hean(q) 


(14) 


Thus  as  m  -*oo,  the  probability  density  function  of  power  scaling  qfc  narrows 

* 

about  Its  mean  value,  while  If  m  -*  0,  It  develops  a  large  cusp  at  the 
origin  and  a  significant  spread  about  the  mean. 


The  normalized  covariance  coefficient  appearing  In  (10)  Is  that  of  the 
power  scaling  variates: 


pw  *  4“  (dk  -  <D(q*  -  q)  •  (15) 

°q 

It  measures  the  degree  of  correlation  between  the  fading  imposed  on  the  signal 
pulses,  l.e.,  frk]  In  (7). 


•This  corresponds  to  a  nonfading  medium. 

10 


The  quantity  K  in  (10)  can  be  Interpreted  as  an  effective  number  of 

t? 

Independent  scalings  in  (4)  and  (7).  For  example,  If  p^  ■  0  for  k  +X, 
then  Kg  *  K;  while  If  *  1  for  all  k,i  ,  then  Kg  =  1.  Both  of  these 
situations  agree  with  physical  Intuition.  No  specific  time  separations  {tkl , 
In  (4)  and  (7),  need  be  assumed  for  (10)  to  apply;  some  of  the  pulses  can  be 
close  together,  while  others  can  be  widely  separated,  subject  of  course  to 
limitation  (3).  The  only  way  that  the  statistical  dependencies  of  the  power 
scalings  {qk]  enter  the  characteristic  function  of  y  is  through  the  double 
summation  of  covariance  coefficients  In  (10). 

However,  It  must  be  remarked  that  the  result  In  (8)  Is  only  an 
approximation,  developed  In  appendix  A.  An  exact  analysis  of  the 
characteristic  function  of  output  y  Is  the  subject  of  a  future  study. 

Exceedance  Distribution  Function  of  Decision  Variable 

Three  alternative  forms  for  the  exceedance  distribution  function  of  y. 


0y(u)  *  Prob(Y  >  u)  ,  (16) 

are  developed  In  appendix  B,  namely,  (B-9),  (B-ll),  and  (B-13),  each  having 
different  merits,  as  discussed  there.  The  one  we  have  used  for  our  numerical 
calculations  Is  the  last  one,  and  Is,  In  fact,  the  detection  probability  when 
signal  Is  present: 
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Is  a  normalized  threshold,  and 


ej<x> 


*"/n! 


Is  the  partial  exponential  [1;  6.5.11].  The  other  two  parameters  appearing  in 
Pp,  namely,  R  and  N,  have  already  been  explained  in  (9)  and  (10). 


The  false  alarm  probability  corresponding  to  (17)  is  obtained  by  setting 


R  =  0: 


Pp  -  exp ( -JO  e^CA) 


Thus,  the  performance  of  the  diversity  combining  processor  Is  governed  by  the 
fundamental  pair  of  equations,  (17)  and  (20).  For  a  specified  number  of 
pulses  K  and  false  alarm  probability  Pp,  (20)  can  be  solved  for  thresholdA. 
These  values  of  K  and  A  are  then  employed  in  (17)  to  evaluate  detection 
probability  PQ  for  any  specified  N  and  E^/N^  A  program  for  this 
procedure  Is  given  In  appendix  C. 
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Detection  Probability  for  Nonfadlnq  Medium 


For  comparison  purposes,  the  detection  probability  for  the  processor  of 
figure  1  In  the  presence  of  a  deterministic,  l.e.,  nonfading,  medium  Is  also 
presented.  This  result  Is  most  easily  obtained  from  the  conditional 
characteristic  function  for  the  system  output,  as  presented  In  (A-13)  and 
(A-14): 

-  (l  -  1  Sufi*  sj.  (21, 

where 


K 

s  -  2.  ak  (22, 

k-1 


Is  the  sum  of  the  power  scalings  on  all  K  pulses.  Since  there  Is  no  fading, 
(21)  Is  directly  the  characteristic  function  of  the  decision  variable. 


A  more  convenient  form  for  (21)  Is  available  when  we  observe  from  (4)  and 
(1)  that  the  received  signal  energy  on  the  k-th  pulse  Is  the  nonrandom  quantity 

2  2  1  12 

Ek  “  A  rk  2  L  *  2  A  L  qk  *  (23) 

Thus,  the  total  received  signal  energy  Is,  using  (22), 

K 

ET  *  2  Ek  "  2  a2l  S  '  (24) 

k-1 
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Coupled  with  the  expression  for  o*  given  in  (9),  we  obtain  for  (21),  the 
characteristic  function. 


f*C)(f) 


In  terms  of  the  parameter 


(;  - 


1$2< 


1  -  1|2oJ  No 


(25) 


_T  =  total  received  signal  energy  over  K  pulses _ 

N  =  single-sided  received  noise  spectral  density  level 


(26) 


The  exceedance  distribution  function  corresponding  to  (25)  is  available 
from  [2;  (50)  and  (51)]  as 


P0-°K 


(27) 


where  we  used  (18)  and  the  Q^-function  defined  as 


QhU.D>  ■[«.  (I)"'1  exp(-  !„.,(«) 


(28) 


The  false  alarm  probability  corresponding  to  (27)  is  obtained  by  setting 
Et  *  0,  and  Is  Identically  (20),  as  it  must  be,  since  the  background  noise 
is  Independent  of  any  signal  fading  characteristics. 


When  we  plot  these  latter  results  for  a  nonfading  medium,  and  compare 

them  with  the  earlier  results  in  (17)  for  the  fading  medium  of  Interest,  we 

replace  ET  by  K  Tj.  This  arbitrary  but  reasonable  assignment  is  necessary 

in  order  to  superpose  the  two  types  of  results  on  one  plot.  However,  strictly 

speaking,  the  parameter  e7/N  appearing  in  the  following  plots  should  be 

1  o 

Interpreted  as  Et/(KNq)  in  the  case  of  no  fading. 
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NUMERICAL  RESULTS 


The  probability  density  functions  (PDF)  for  five  values  of  the  scaling 
parameter  m  in  (11)  and  (13),  representing  five  distinct  possible  fading 
behaviors,  are  plotted  In  figure  2.  All  densities  of  the  class  considered 
vanish  for  negative  argument.  So  these  behaviors  are,  for  the  amplitude 
variate: (1)  an  Impulse  representing  the  nonfading  case,  which  Is  also  the 
limiting  case  of  the  class  as  m  -»o6;  (2)  continuous  with  continuous  derivative 
at  the  origin  (m  =  2);  (3)  continuous  with  discontinuous  derivative  at  the 
origin  (m  =  1 ,  Rayleigh  amplitude  density);  (4)  a  finite  discontinuity  at  the 
origin  (m  =  1/2,  a  one-sided  or  full-wave  rectified  zero-mean  Gaussian 
density);  and  (5)  an  Infinite  cusp  at  the  origin  (m  =  1/4).  The  amplitude 
scaling  r,  with  the  probability  density  function  of  (11),  has  second  moment 

r2  *  mrf . 

The  normalized  amplitude  scaling 

*r-r /(?Yn. 

normalized  by  the  root-mean-square  (rms)  value  of  the  amplitude  scaling  (not 
the  standard  deviation  or),  has  probability  density  function 

p?(u)  =  -2-m  u  p(m)Xpt~mu2*  for  u  >  0  .  (31 

For  the  five  values  of  m  (  =  *,  2,  1,  1/2,  1/4),  the  probability  density 
functions  for  this  rms -norma  1 i zed  amplitude  scaling  are  shown  in  figure  2A. 


(29) 


(30) 
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The  corresponding  probability  density  functions  for  power  scaling  q  are 
chi-squared  variates  with  2m  degrees  of  freedom,  as  given  in  (13).  The  mean 
of  variate  q  is 


q  *  m* . 

The  normalized  power  scaling 


(32) 


q  *  q/q  ,  (33) 

normalized  by  the  mean  power,  has  probability  density  function 

PS(U)  =  ^  for  u  >  0  .  (34) 

The  normalized  power  scaling  probability  density  functions  are  plotted  in 
figure  2B  for  the  same  five  values  of  m. 


For  each  finite  value  of  m  above,  and  for  each  of  the  four  combinations 
of  system  performance  parameters  Pp  =  .5  or  .9,  and  Pp  =  IE-6  or  IE-8, 
the  required  minimum  system  input  signal-to-noise  ratio  (SNR)  to  achieve  this 
performance,  as  a  function  of  the  number  of  pulses  K  in  a  pulse  train,  has 
been  computed  using  the  computer  programs  listed  in  appendix  C  which  implement 
the  theoretical  analysis  above.  Although  this  analysis  and  these  programs  can 
handle  general  pulse-to-pulse  normalized  covariance  coefficients  {p^}  .  The 
figures  have  been  computed  for  the  special  case  of  exponential  correlation 

Ik-il 


p^  =  p 


-1  <  p  <  1  . 


(35) 
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Each  of  the  graphs  has  two  sets  of  curves:  the  solid  curves  are  for  the 
required  average  received  signal  energy  per  pulse  (T!j /Nq)  in  dB,  as  a 
function  of  the  number  of  pulses  K,  and  the  dotted  curves  are  for  the  total 
average  received  signal  energy  (E^./N0)  in  dB.  For  ease  of  comparison, 
each  set  of  curves  on  each  graph  includes  its  nonfading  (NF)  or  m  -m» 
counterpart.  In  every  case  computed,  the  solid  curves  show  that  the  energy 
per  pulse,  or  input  signal-to-noise  ratio  required  for  acceptable  system 
performance,  decreases  as  the  number  of  pulses  increases.  This  matches  the 
intuitive  argument  that  using  more  pulses  implies  we  are  using  more  total 
energy,  so  we  should  be  able  to  maintain  the  same  level  of  performance  with 
less  energy  per  pulse.  The  Intuitive  conclusion  is  correct,  but  scanning  the 
dotted  sets  of  curves  (representing  total  energy  in  the  pulse  train)  shows 
this  intuitive  argument  is  only  part  of  the  story.  For  more-or-less 
uncorrelated  adjacent  pulse  fading  and  for  the  smaller  scaling  parameters 
(m  <  1),  the  total  Input  signal  energy  required  initially  decreases  as  we 
Increase  the  number  of  pulses  In  a  pulse  train.  For  a  large  number  of  pulses, 
the  total  energy  will  eventually  increase  with  K,  but  for  a  small  number  of 
pulses,  a  short  burst  may  be  more  efficient.  For  example,  with  =  .5, 

Pp  =  IE-8,  and  a  one-sided  Gaussian  amplitude  probability  density  function, 
m  =  1/2,  figure  13  shows  that  the  total  energy  required  reaches  a  minimum  of 
14.4  dB  for  p  =  0  with  3  pulses  per  burst;  for  p  =  .5,  the  minimum  of  15.2  dB 
occurs  with  5  pulses.  In  all  cases,  the  p  =  .5  curve  is  closer  to  the  p  =  0 
curve  than  to  the  p  =  1  curve,  and  the  high  correlation  p  =  1  curve  always 
Increases.  Thus,  for  moderate  values  of  correlation,  we  can  expect  system 
performance  to  be  robust  and  system  performance  will  be  severely  degraded  only 
for  highly  correlated  fading. 
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As  the  exponential  covariance  coefficient  p  in  (35)  decreases  toward  -1, 
the  performance  curves  for  even  K  approach  those  corresponding  to  the 
nonfading  case.  This  is  easily  seen  since  the  sum  of  two  successive  power 
scalings  In  the  sum  S  In  (22)  and  (A-14)  then  tends  to  2q,  twice  the  mean  of 
q,  which  Is  a  nonrandom  quantity.  Equivalently,  we  can  say  the  random 
fluctuations  In  the  power  scalings  {qj  have  canceled  each  other.  For  K 
odd,  however,  there  Is  always  one  unpaired  random  fluctuation,  even  for  P  =  -1 ; 
this  means  that  for  K  odd  and  p  *  -1,  system  performance  is  slightly  poorer 
than  for  the  nonfading  situation.  As  K  Increases  through  odd  values,  this 
discrepancy  decreases  to  Insignificance.  Even  for  K  *  3,  the  difference  is 
small.  This  has  been  confirmed  analytically  and  verified  by  numerical 
calculation  for  values  of  p  very  near  -1,  and  for  several  K;  these  results 
have  been  omitted  from  the  figures  which  include  only  nonnegative  values  of  p. 

A  separate  Issue  Is  whether  values  of  p  near  p  -  -1  can  be  realized. 

That  Is,  if  the  first-order  probability  density  function  of  the  nonnegative 
power  scaling  variate  q  is  not  symmetric  about  Its  mean  q,  then  arbitrarily 
large  negative  values  of  p  are  not  allowed.  For  example,  for  p^  to  reach 
-1 ,  we  would  need  (qfc  -  q)  »  -(q^  -  q)  with  probability  one.  However,  for 
the  class  of  power  scalings  treated  here,  this  is  clearly  impossible:  large 
positive  values  of  q^  -  q  occur  with  nonzero  probability,  for  which  the 
corresponding  negative  values  -(qfc  -  q)  occur  with  probability  zero.  Thus, 
for  most  probability  density  functions,  and  In  particular  for  the  class 
assumed  here,  the  p  =  -1  results  would  represent  an  unreachable  bound  to 
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system  performance.  This  Is  an  Illustration  of  the  fact  that  the  first-order 
probability  density  function  and  the  correlation  (or  spectrum)  of  a  random 
process  cannot  be  Independently  specified  in  general.  Complete  first-order 
Information  Imposes  some  restrictions  on  allowable  second-order  statistics. 
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Input  SNR  in  dB  Input  SNR  in  dB 
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SUMMARY 

One  manner  of  summarizing  the  behavior  of  the  amplitude  and  power  scaling 
variates  (figure  2)  Is  to  observe  that  as  the  parameter  m  decreases  from 
Infinity  (the  nonfading  case),  the  fading  becomes  deeper  more  often.  For 
small  m,  the  very  deep  fades  occur  very  often.  As  m  nears  zero,  the  fades  are 
so  deep  and  occur  so  often  that  large  arguments  of  the  scaling  probability 
density  functions  are  rarely  observed.  Thus,  for  a  given  level  of 
performance,  the  required  Input  slgnal-to-nolse  ratio  must  increase  as  m 
decreases. 

This  conclusion  Is  borne  out  by  the  Pp  *  .9,  Pp  **  1 E— 6  case 
(figures  4,  8,  12,  and  16)  where,  for  a  single  pulse  system,  K=l,  the  required 
Input  slgnal-to-nolse  ratio  Increases  from  17  dB  to  47  dB  as  m  decreases  from 
2  to  .25.  Relaxing  the  probability  of  detection  to  PQ  =  .5  relaxes  the 
corresponding  slgnal-to-nolse  ratios  to  the  range  from  12  dB  to  19  dB.  It  is 
also  not  surprising  that  tightening  the  probability  of  false  alarm  to  IE-8 
further  Increases  the  required  slgnal-to-nolse  ratio.  For  m  =  .25  and 
PF  =  IE-8,  the  slgnal-to-nolse  ratio  requirements  are  preposterous  and  the 
corresponding  figures  for  Pp  =  .5  and  Pp  =  .9,  which  would  have  completed 
the  set,  have  been  omitted. 

Although  the  signal  model  employed  here  In  (1)  and  (4)  utilizes  a  common 
carrier  frequency  f  for  all  K  pulses,  the  results  extend  immediately  to  the 
case  of  arbitrary  different  carrier  frequencies  ffe  for  the  k-th  pulse, 
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provided  that  these  various  frequencies  are  known  and  appropriately  used  at 
the  receiver.  Also,  although  the  particular  signal  pulse  In  (1)  is  a  simple 
burst  of  a  slnewave,  the  current  results  extend  to  arbitrary  signal  waveforms 
Pk(t)  on  each  pulse,  provided  that  matched  filters  are  employed  in  the 
receiver,  and  that  each  waveform  ( t)  have  identical  transmitted  energy. 
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APPENDIX  A 

CHARACTERISTIC  FUNCTION  OF  SYSTEM  OUTPUT 


Hatched  Filter  Output  Characterl zatlon 


In  this  appendix,  we  compute  the  characteristic  function  of  the  decision 
variable,  described  in  (7)  et  seq.  The  In-phase  and  quadrature  outputs  of  the 
k-th  matched  filter,  at  the  time  instant  of  peak  signal  output,  are 


j^k  *  dt  |a  rk  cos(2*fQ(t  -  tk  -  D)  +  ©k)  +•  n(t)j 


*  [sin](2,rfo(t  “  tk  "  D))  for  1  <  k  <  K  ,  (A-l) 

where  use  has  been  made  of  (1)  and  (7),  and  phase  shift  eQ  has  been  absorbed 
In  random  phase  shift  ©k,  with  no  loss  of  generality.  That  is, 

Ik  =  <*k  +  **k  *  Qfc  =  ®k  +  wk  ’ 

where  the  signal  components  are,  upon  use  of  the  facts  that  If  »  1  and 

o 

that  there  is  no  overlap,  as  indicated  in  (3), 


*k  *  2AL  rk  cos  8k  . 


°k  ’  -  2Al  rk  s1n  ®k  • 


and  the  noise  components  are 


n<t» 


2*f0(t  -  tk  -  D» 


( A— 3 ) 


( A— 4 ) 
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The  noise  components  in  (A-4)  are  Gaussian  zero-mean  random  variables  since 


n(t)  =  0  .  n(t1)n(t2)  =  «(t1  -  t2) 


( A— 5 ) 


There  also  immediately  follows  by  the  use  of  (A— 5) . 


2  2  M  ,  1  „  2  -  n 

*k  =  wk  =  Nd  L  2  2  °n  •  Vk  =  0 


( A— 6 ) 


Furthermore,  due  to  the  lack  of  overlap,  as  in  (3),  the  random  variables 
vk,»k  are  independent  of  for  k  #  j.  Thus,  the  collection  of  2K 

random  variables  in  (A-4)  are  all  independent  and  identically  distributed. 


Conditional  Characteristic  Function 


The  k-th  squared  envelope  sample  at  the  matched  filter  output  is 


\  -  ijj  +  Qk  -  («k  +  wk)2  ♦  (Bk  +  \)2  . 


( A— 7 ) 


by  reference  to  (A-l)  and  ( A— 2 ) .  We  will  temporarily  consider  that  the  signal 
components  in  (A-3)  are  fixed  non-random  constants;  then  the  conditional 
characteristic  function  of  random  variable  is,  uppn  use  of  (A-6), 


fk(?)  =  exp( if Yk)  =  exp[i?(ak  +  pk)2  +  1$ (Dk  +  «k ) 2 ] 


=  (Tdp  dv  ——  exp/-  * — exp [ i"J ( 

JJ  >'°l  ^  2°ll 


<*l,  +  u)2  +  i|(B.  +  v)2] 


( A -8) 


Now  let 


=  p  cos  e  ,  v  =  p  sin  e  , 


(A  -9) 


&IXKS 
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to  get 


-  (1  -  l^ojj)"1 


exp 


W  * 8') 

-  '**4  .  ' 


Here  we  employed  [3;  8.411  1  and  6.631  4]. 


( A— 1 0) 


The  output  of  the  system  Is  the  sum  of  the  squared  envelope  samples  In 


( A— 7 ) : 


Y 


<‘k  *  Qk>  • 


(A-l 1 ) 


Since  all  the  noise  random  variables  in  (A-4)  are  Independent,  the  conditional 
characteristic  function  of  system  output  random  variable  y  follows  from  ( A-l 0) 
as 


f<C)<?)  «  (1  -  1f2oj;)'K  exp 


if 


1  -  1?2o' 


<«k  +  Bk> 


n  k=l 


(A-l  2) 


If  we  now  use  (A-3),  this  can  be  expressed  as 


^°(?)  -  (1  -  1*2oj;)"K  exp 


-Li —  m 


1  -  If 2. 


2  V  2/ 


2 

k=l 


( A— 1 3 ) 
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It  Is  important  to  observe  that  this  conditional  characteristic  function 
depends  on  the  signal  amplitude  scalings  ^r^  in  (4)  and  (7)  only  through 
the  sum  of  squares 


‘.V»1 


f 

i 


where 


qk  =  rfe  *  power  scaling  of  k-th  pulse 


First  Order  Statistics  of  Scallnos 
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This  Is  recognized  as  the  probability  density  function  of  a  chi-squared  random 
variate  with  2m  degrees  of  freedom;  however,  2m  need  not  be  an  Integer  here. 
The  general  »-th  moment  of  q^  Is 


v  v  r*(m  +  w) 
qk  “  a  T(m) 


and.  In  particular,  we  have 


( A— 1 9) 


r  _  2  Std  Dev(a)  1 

q  *  am,  Var(q)  =  a  m,  - — — 

q  Vm 


(A-20) 


The  first  order  characteristic  function  of  power  scaling  qfc  is  given  by 


-m 


the  Fourier  transform  of  (A-18): 

fq(f)  -  j*  du  exP(1?u)  pq(u)  »  (1  -  if«)' 

First  Order  Characteristic  Function  of  Sum  S  In  Special  Cases 


( A— 21 ) 


If  all  the  power  scalings  {qfc]  in  sum  S  In  (A-14)  were  Independent, 
then  the  characteristic  function  of  S  would  be  given  by 

(1  -  1X«)"mK  ,  ( A— 22 ) 

as  may  be  seen  from  (A-21).  This  situation  could  come  about  approximately 
when  the  signal  pulses  In  (2),  (4),  and  (7)  are  widely  separated  in  time  and 
subject  to  uncorrelated  (linearly  Independent)  fading. 
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On  the  other  hand,  If  all  the  power  scalings  were  completely  dependent. 


then  the  characteristic  function  of  sum  S  would  be 

(1  -  IfaK)""’  ( A— 23) 

since  S  =  Kq1  In  this  case;  see  ( A— 1 4)  and  ( A— 21 ) .  This  case  corresponds  to 
close  signal  pulses  and/or  very  slow  fading. 

Approximate  First  Order  Characteristic  Function  of  Sum  S 

* 

The  latter  two  results  suggest  the  following  form  as  an  approximation 
to  the  characteristic  function  of  sum  S  In  (A— 14)  when  the  power  scalings 
[qj  are  partially  dependent: 

-mK 

fs(f)  *  (1  -  1?B)  e  .  ( A -24 ) 

This  latter  form  of  characteristic  function  has  mean  mK  B  and  variance 

e 

2 

mKgB  .  We  will  choose  the  two  unknown  parameters  in  ( A— 24 )  so  that  these 
two  statistics  Identically  equal  the  corresponding  exact  values  determined 
directly  from  (A-14).  Specifically,  using  (A-20),  mean 

S=Kq=Kam,  (A-25) 


*A  similar  procedure  was  employed  with  great  success  for  a  spectral  analysis 
technique  In  [4;  (38)  et  seq.]. 
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and  variance 


where 


Var(S)  -  (S  -  S)2 


K  K 

°q  -  <*2">  . 

k  Jl  =1  kj=l 


( A— 26 ) 


“  “2  (qk  '  q)(qi  “  q)  (A-27) 

°q 

Is  the  normalized  covariance  coefficient  of  power  scalings  fqj.  Combining 
these  quantities,  we  find  that 


2L  pkji 

kjM 

6  -  -  r  ■  S  r  • 

e  e 

The  quantity  Kg  In  (A-24)  and  (A-28)  can  be  interpreted  as  an  effective 
number  of  Independent  scalings  in  (A-14).  For  example.  If  p ^  =  0  for 
k  *Ji,  then  Kg  =  K;  while  If  p^  =  1  for  all  k,^,  then  Kg  =  1.  Both  of 
these  situations  agree  with  physical  Intuition.  Also,  from  (A-28),  B  =  a  in 
the  former  case,  and  Q  =  aK  In  the  latter  case;  (A-24)  then  reduces  to  (A-22) 
and  (A-23),  respectively,  as  required. 
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No  specific  time  separations  ft^},  In  (4)  and  (7),  need  be  assumed  for 
(A^28)  to  apply.  Some  of  the  pulses  can  be  close  together,  while  others  can 
be  widely  separated,  subject  of  course  to  limitation  (3).  The  only  way  that 
the  statistical  dependencies  of  the  power  scalings  {q^  In  sum  S  enter  the 
characteristic  function  of  S  Is  through  the  double  summation  of  covariance 
coefficients  In  (A-28). 

Unconditional  Characteristic  Function  of  Output  v 


The  first  order  probability  density  function  corresponding  to 
characteristic  function  (A— 24)  is 


Ps(u) 


T(N)  0 


N 


for  u  >  0 


(A -29) 


where 


N  =  m  Ke  .  (A-30) 

None  of  the  parameters,  N,  m,  Kg,  need  be  Integer.  Also,  N  can  be  larger  or 
smaller  than  K,  the  number  of  signal  pulses. 


The  conditional  characteristic  function  of  system  output  y  was  given  in 
(A-13).  The  unconditional  characteristic  function  of  y  Is  obtained  by 
averaging  (A-13)  with  respect  to  (A-29): 


38 


TR  7707 


oo 

f  (f)  -  f  du  ^  (1  - 

«  J  ft" 


r(N)  8 


1J2 a2)~K  exp 


fr^ : 


(1  -  it2a2n)H~K  (1  -  1J2<4)~N  , 


(A-31) 


where 


To  obtain  this  last  relation,  we  used  (A-6),  (A— 28) ,  (A-30),  (6),  and 

Nq  -  2Nd  .  ( A-33) 

the  single-sided  noise  density  level  stated  just  above  (7). 


Equation  (A-31)  Is  a  compact  simple  expression  for  the  characteristic 
function  of  decision  variable  y.  It  could  be  used  directly  In  the  efficient 
procedures  of  [5]  to  get  the  cumulative  and/or  exceedance  distributions  of  y. 
In  the  next  appendix,  we  derive  analytic  expressions  for  these  probabilities. 
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APPENDIX  B 


EXCEEDANCE  DISTRIBUTION  OF  SYSTEM  OUTPUT 


The  characteristic  function  of  decision  variable  y  is  given  in  (A-31). 


In  order  to  find  the  corresponding  exceedance  distribution  function,  we  use 


the  procedure  in  [6;  appendix  A]:  since  we  know  the  characteristic  function- 


probability  density  function  pair 


(1  -  iya)J 


uJ  1  exp(-u/al 


T(J)  aJ 


for  u  >  0,  (J  >  0)  .  ( B— 1 ) 


the  probability  density  function  corresponding  to  multiplicative 


characteristic  function  (A-31)  is  given  by  convolution 


U 

»T(U)  . 


xJ-1  expf-x/al 


:u  -  x] 


T(J)  a‘ 


»xp(-(u  -  x)/t 
l/U\  L.  N 


P(N)  b 


for  u  >  0  ,  (B-2) 


where 


J  -  K  -  N.  a  -  2^.  b  -  2oj 


( B— 3 ) 


We  presume  N  <  K  for  now,  in  order  that  J  >  0. 


Employing  [3;  3.383  1  and  8.384  1],  (B-2)  becomes 


P  (u)  -  U 

t  a*""  bN  f(K) 


lFl(^:  K;  u(a  "  b))  for  u  >  °-  (B"4) 


Since  (A-31)  is  an  analytic  function  of  N,  as  is  its  transform  (B-4),  we  can 


now  analytically  continue  (8-4)  to  values  of  N  >  K, 
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The  exceedance  distribution  function  of  y  Is 


V, 

g 


Prob(Y  >  u)  -  0^(u) 


fix  pY 


(x)  = 


y 

,y 


P(K)  (b/a)' 


oo 

J  “  «_t  t"'1  S-1  0 


( B— 5) 


I 


4 

A. 

*>jj 

|VJ 


I 


Reference  to  (8-3)  and  (A-32)  yields 


Qv(u)  -  - 

y  r(K) 


00 

77777]“  *rht)- 


where 


D  1  K  .  u 

R  *  N  N  *  y  2  • 

0  2a 


(8-6) 


(B-7) 


By  expanding  ^  In  a  power  series  according  to  [3;  9.210  1],  Integrating 
term  by  term,  and  using  the  partial  exponential  [1;  6.5.11] 


ej(x)  *  ^  xn/n!  , 


(8-0) 


(8-6)  develops  Into  the  form 


v“> 5^  W 


Vi«w  • 


(8-9) 


:.a 


This  series  for  the  exceedance  distribution  function  converges  for  all  R,  but 
rather  slowly  for  large  R. 
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For  signal  present,  (B-9)  gives  the  detection  probability.  When  signal 
Is  absent,  R  =  0  from  (B-7),  and  (B-9)  reduces  to  the  false  alarm  probability 

Pp  *  exp(-A)  e^_-|  (-*■)  ,  (B-10) 

which  Is  Independent  of  N  (defined  In  (A-30)  and  (A-28)). 


An  alternative  form  to  (B-9)  Is  available  If  we  employ  Rummer's 
transformation  [1;  13.1.27]  In  (B-6): 


OO 


0Y(u) 


<’ +  R>K'N  1  (-»>"  eK-u„(rh)-  <»-"> 

n=0 


If  N  Is  an  Integer  and  If  N  >  K,  this  series  terminates,  and  Is  composed  of 
all  positive  terms,  which  makes  It  very  attractive.  However,  if  either 
condition  Is  violated,  then  (B-11)  converges  only  for  R  <  1,  but  diverges  for 
R  >  1.  Since  this  latter  case  Is  of  practical  significance,  (B-11)  is  then 
not  too  useful  for  general  values  of  N. 


The  only  convergence  factor  in  (B-9)  is  the  R/(l  +  R)  term.  However,  we 
can  create  another  convergence  term  by  adding  and  subtracting  exp(A)  from  the 
partial  exponential,  since  this  Is  Its  limit  value  as  n  tends  to  Infinity;  see 
(B-8) .  By  then  employing  the  result  that  [3;  9.100  and  9.121  1] 

j>  8:  rr-p)-  C  ♦  *>"  •  o-12* 

n=0 
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(8-9)  becomes 


V“>  - 1  -  ^  (rh)n  f*«pw  -  v,«w]  • 


(B-13) 


This  last  form  was  used  for  all  the  numerical  results  obtained  here;  a  program 
for  this  expansion  is  given  In  appendix  C.  If  necessary,  a  more  elaborate  and 
accurate  procedure  Involving  only  positive  summations  Is  developed  in 
[6;  pages  A-6  to  A-8]. 


To  summarize,  the  parameters  required  in  detection  probability  (B-13)  are 


-A  -  — jj-  .  «  *  S1  N  *  (B-14) 

2  aj,  0 

n 

2 

where  an  is  given  by  (A-6)  and  N  is  given  by  (A-30)  and  (A-28).  Since  A 
Is  simply  a  scaled  version  of  u,  we  can  equally  well  keep  JV  as  the  fundamental 
threshold  variable  In  (8-10)  and  (B-13). 
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PROGRAM  LISTINGS 

Two  programs  are  listed  here,  the  first  for  fading,  the  second  for 
no-fading,  of  the  received  signal.  Inputs  required  of  the  user  for  the  fading 
case  are  given  in  lines  20  -  50: 

Pf  probability  of  false  alarm 

Pd  probability  of  detection 

M  parameter  m  in  (11)  and  (13) 

Rho  =  p'*-*'  in  (10)  and  (15)  . 

The  particular  exponential  normalized  covariance  p^  programmed  in  line  1490 
can  be  easily  replaced  by  other  more  general  cases  of  interest  to  the  user. 

The  output  of  the  main  program,  line  290,  Is  the  value  in  dB  of  FJ/Nq 

In  (9),  namely  a  measure  of  the  required  system  Input  slgnal-to-noise  ratio  to 
achieve  the  desired  false  alarm  and  detection  probabilities.  The  number  of 
pulses,  K,  Is  taken  to  be  1  to  10  In  line  120,  but  can  be  easily  changed.  The 
first  Inverse  function  subroutine  in  line  770  solves  for  the  required 
threshold  for  specified  Pp  and  K.  The  second  subroutine  in  line  1130  solves 
for  the  required  1T.|/No  for  specified  PQ  and  N. 

A  similar  explanation  holds  for  the  no-fading  program  listed  at  the  end 

of  the  appendix.  Duplicate  routines  are  not  listed  In  full,  in  order  to  save 

space.  The  word  DOUBLE  denotes  INTEGER  variables  in  Hewlett  Packard  BASIC  on 
the  9000  computer. 
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10  '  3  N  R  -  F  A  D I N  G 

20  Pf = 1 . E-6 

30  Pd=. 5 

40  M=l. 

50  Rho=.5 

60  A*=" IE-6, .5,  1 .  ,  .5" 

70  Ef=l.E-20 

30  Ed=  1 .  E-  1 5 

90  DIM  ElnodbC 1 ; 10> 

100  COM  DOUBLE  K,REAL  N,Thr 

110  DOUBLE  Ks,  Ls 

120  FOP  K  =  1  TO  10 

130  3=0. 

140  FOP  K s  = 1  TO  K 

150  FOR  L s  = 1  TO  K 

160  3  =  3  +  F  N  C  o  n  o  r  m  i  1 

170  NEXT  Ls 

130  NEXT  Ks 

190  Keff=K*K-S 

200  N»M*Keff 

210  Thrineal. 

20  T  h  r  s  t.  ir  t  =  - 

_30  CALL  I  n  •>  e  r 

240  E  1  ri  o  3 1-  ar  t  =  1 0 

250  Elnoinc=l. 

260  CftLL  Inc 

270  ElnodbL  K 

230  NEXT  K 

290  PRINT  E 1 nodb < * > 

300  CREATE  DATA  A*,  4 

310  ASSIGN  #1  TO  A* 

320  PRINT  #1 ; Elnodb' *> 

330  ASSIGN  #1  TO  * 

340  END 

350  ' 

360  DEF  FNPf  < Thr 

370  COM  DOUBLE  K 

330  DOUBLE  Ks 

390  3=T=EXP< -Thr > 

400  FOR  K 3  =  1  TO  K -  1 

410  T=T*Thr  Is 

420  3  =  3  ♦  T 

4  30  NEXT  K  s 

440  RETURN  -3 

450  FNEND 

4  6  0  1 


FALSE  ALARM  PROBABILITY 

DETECTION  PROBABILITY 

FADING  PARAMETER  < SINGLE  PULSE  > 

FADING  CORRELATION  COEFFICIENT 

Pf ,Pd,M, Rho 

TOLERANCE  ON  Pf 

TOLERANCE  ON  Pd 

REQUIRED  INPUT  SNR  PER  PULSE  < D B ) 


NUMBER  OF  PULSES 


THRESHOLD  INCREMENT 
!  THRESHOLD  STARTING  VALUE 


,  L  3 ,  P  t'l  o 


L  0  G  <  P  f  :>  -  T  h  r  i  n  c 

3  f  u  n  c  *,  i  o  n  1  (  -  P  t  ,  £  f  ,  T  h  r  s  t ar  t  ,  T  h  r  l  n  c  ,  T  h  r  » 

!  El  No  STARTING  VALUE 
'  El -No  INCREMENT 

•  r  sf  uric  t.  i  on2 <  Pd ,  Ed ,  E 1  nos t  art. ,  E 1  no  i  nc  ,  E 1  nc 
■  =  10. *LGT  < E 1  no  > 


FALSE  ALARM  PROBABILITY 


'  -  TO  YIELD  INCREASING  FUNCTION 
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Cop? 


to  Dlic  does  not 

iuliy  ieflible  repioductioa 


QV 


il  M 


DEF  FNPd  <  E 1  no >  !  FADING  DETECTION  PROBAB I  LIT’ 

COM  DOUBLE  K , REAL  N,Thr 

Error® 1 .  E- 10 

DOUBLE  K 1 ,  Ks 

R=Elno*K/N 

Et=EXP<Thr> 

K 1 =  K- 1 
N1=N-1. 

R 1 = 1 . +R 
Q=R  R1 
E®Te=  1 . 

FOR  Ks®  1  TO  K 1 
Te=T**Thr  /Ks 
E=E+Te 
NEXT  Ks 

S»B»MAX<Et-E,0.  > 

T=  1  . 

FOR  Ks® 1  TO  1000 
Te®T e*Thr  ' < K 1 +  Ks > 

B=MAX<B-Te, 0. > 

T  =  T  *  Q  *  <  N 1  +  k  s  >  /  K  s 

P  r  =  T  *  B 

3=S*Pr 

IF  AB3  <  Pr  >  <  =Err or*ABS  (•$  >  THEN  730 
NEXT  Ks 

PRINT  "1000  TERMS  AT: " ; Kj N; Thr 5 Elno;  Pr-  S 
P  d  =  1 . -EXP  -Thr-N*L0G<R1 > >*S 
RETURN  Pd 
FNEND 

I 

SUB  I  n v erst' unc  t  i  on  1  <  De  s  i  red ,  Error ,  X 1 ,  Be 1  ,  X2  > 

X2*X 1+D* 1 
F  1  ®FNPf  <  X 1 ') 

F2  =  FNPf  <  X2  > 

IF  F2>®Desired  THEN  8e© 

X 1  =  X  2 
X2=X2+De 1 
F  1  =F2 
GOTO  800 

IF  Fl<*Desir*d  THEN  92© 

X  2  ®  X  1 
XI *x 1 -De  I 
F2  =  F  t 

Fl=FNPf C XI  > 

GOTO  S3© 

X  »=  X  1 
Xb=X2 

IF  F2-Des i red : Des i red-Fl  THEN  101© 

T  «  X 1 

x  i  =:  :2 

x  2 = T 

T=F  1 
F  1  =  F  2 
F2  =  T 


Copy  available  to  DTIC  does  not 
penaii  fully  legible  reproduction 


1010  IF  RBSCF2-Desired><Error  THEM  1110 

1020  IF  F2  =  F  1  THEN  1110 

1030  T*<X1*<F2-Dest  red>-X2*,:Fl-D**i  red  >  >  s  <  F2-F 1  > 
1040  T=MRX<T,Xa> 

1050  T*M I N <  T , Xb  > 

1060  X 1 *X2 

1070  X2  =  T 

1080  F 1 =F2 

1080  F2=FNPf<X2> 

1100  GOTO  1010 

1110  SUBEND 

1120  ! 

1130  SUB  Inoersfunc t  i  ort2<Des i  red,  Error ,  XI ,  De  1  ,  X2> 
1140  X2*Xl+De 1 

1150  Fl»FNPd<Xl> 

1160  F2=FNPd<X2> 

1170  IF  F2>=Des i red  THEN  1220 

1180  X1-X2 

1190  X2*X2+De 1 

1200  F 1 =F2 

1210  GOTO  1160 

1220  IF  F 1 <  *Des l red  THEN  1280 

1230  X2-X1 

1240  X 1 *X 1 -De 1 

1250  F2=F1 

1260  F 1 =FNPd<  X 1 > 

1270  GOTO  1220 

1280  Xa=Xl 

1290  Xb«X2 

1300  IF  F2-Desi redCDesi red-Fl  THEN  1370 

1310  T =X 1 

1320  X 1 »X2 

1330  X2*T 

1340  T=F 1 

1350  F 1 =F2 

1360  F2=T 

1370  IF  PBS <F2-De*i red >< Error  THEN  1470 

1380  IF  F2=F1  THEN  1470 

1 390  T«<X1*<F2-De#i  red) -X2*< F 1  -De* i  red  >  >  •••" <F2-Fl  > 

1400  T = M ft X  «r  T,  Xa) 

1410  T*M I N<  T , Xb  > 

1420  X 1 =  X  2 

1430  X2*T 

1440  F 1 *F2 

1450  F2  =  FNPd  < X2 ) 

1460  GOTO  1370 

1470  SUBEND 

1480  ! 

1490  DEF  FNCoMriorni 1  DOUBLE  K  £ , Ls , RERL  Rho  • 

1500  RETURN  Rho  RBS- Ks-L*  ' 

1510  FNEND 


Copy  avnil  ’  !  •  to  DTIC  does  not 

penmt  fully  legible  reproduction 
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10  '  SNR-NO-FAD I NG 

20  Pf *  1 .  E-6  ! 

30  Pd=.3  ! 

40  A*=" IE-6, . 5"  ! 

30  £f*l.E-20  ! 

60  Ed* 1 . E- 1 3  ! 

70  DIM  ElnodbC 1 : 1 0 >  ! 

30  COM  DOUBLE  K , REAL  Thr2 

90  FOR  K* 1  TO  10  ! 

1O0  Thrinc*l.  ! 

110  Thrstart.*-LOG<Pf  >-Thri  nc 

120  CALL  Inver sf unc t i on  1 < -Pf 

130  Thr2=SQP<2.*Thr> 

140  Elnoit art=10.  ! 

150  El  no >nc=l.  * 

160  CALL  Inwerifunct i on2<Pd, 

170  ElnodbC  10*10.  *LGT<Elno> 

130  NEXT  K 

190  PRINT  Elnodb<*> 

200  CREATE  DATA  A*,  4 

210  ASSIGN  #1  TO  A# 

220  PRINT  #l;Elnodb<*> 

230  ASSIGN  #1  TO  * 

240  END 

250  ! 

260  DEF  FNPfCThr)  ! 

270  COM  DOUBLE  K 

230  DOUBLE  Ks 

290  S*T*EXP  < -Thr > 

300  FOR  Ks* 1  TO  K- l 

310  T*T*Thr/K* 

320  S  *  S  +  T 

330  NEXT  Ks 

340  RETURN  -S  ! 

330  FNEHD 

360  1 

370  DEF  FNPd < E 1  no >  • 

380  COM  DOUBLE  K, REAL  Thr 2 

3 9 0  P d * F N Q rn  <  K  ,  S Q R 1  2 .  *K*E  1  n ) 

400  RETURN  Pd 

410  FNEND 

420  1 


FALSE  ALARM  PROBABILITY 
DETECTION  PROBABILITY 
Pf  ,Pd 

TOLERANCE  ON  Pf 
TOLERANCE  ON  Pd 

REQUIRED  INPUT  SNR  PER  PULSE  <.DB> 

NUMBER  OF  PULSES 
THRESHOLD  INCREMENT 
!  THRESHOLD  STARTING  VALUE 
Ef , Thrst  art , Thr i nc , Thr  > 

El -'No  STARTING  VALUE 
El -"No  INCREMENT 
d, Elnost  art , Elnoi nc , Elno) 


FALSE  ALARM  PROBABILITY 


-  TO  YIELD  INCREASING  FUNCTION 


NO-FADING  DETECTION  PROBABILITY 
Thr 2  > 


Copy  available  to  DTIC  does  n<fc 
permit  full?  legible  reproduction 
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43©  DEF  FNQrnC  DOUBLE  M  ,  REAL  R  ,  B  >  !  QM  <  A,  B  > 

440  Err or* 1 . E- 1 7 

430  DOUBLE  M1,J 

460  Q3*.5*A*A 

470  Q4= . 5*B*B  ' 

480  Q5=EXP<-.5*<Q3*Q4>:> 

49©  06*07*05 

500  M 1 =  M- 1 

510  FOR  J*1  TO  Ml 

520  Q7=07*Q4-"J 

530  06=06+07 

540  NEXT  J 

550  0m=Q5*Q6 

560  FOR  J=1  TO  1000 

570  05=05*03- J 

530  07=07*04'-' <  J  +  Ml  > 

590  06=06*07 

600  09=05*06 

610  0m=0m*Q9 

620  IF  Q9<  =Er  ror*Qni  THEN  650 

630  NEXT  J 

640  PRINT  "1000  TERMS  IN  FNQm < N,  A. B  -  AT  M;M;A?B 

650  RETURN  Qm 

660  FNEND 

670  ' 

SUB  Irwersf  unct  i  oril  (  Desi  red,  Error ,  XI ,  Del  ,  X2 

'  LISTED  ABOVE 

SUB  I r,uersf  unc t i on2 1 Des i red, Error , X 1 , De 1 , X2 

'  LISTED  ABOVE 
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